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The dispersion relations for plasmons in intrinsic and extrinsic semimetallic armchair graphene 
nanoribbons (acGNR) are calculated in the random phase approximation using the orthogonal 
Pz-orbital tight binding method. Our model predicts new plasmons for acGNR of odd atomic 
widths A^ = 5,ll,17, ... Our model further predicts plasmons in acGNR of even atomic width N = 
2, 8, 14, . . . related to those found using a Dirac continuum model, but with different quantitative 
dispersion characteristics. We find that the dispersion of all plasmons in semimetallic acGNR 
depends strongly on the localization of the pz electronic wavef unctions. We also find that overlap 
integrals for acGNR behave in a more complex way than predicted by the Dirac continuum model, 
suggesting that these plasmons will experience a small damping for all ^ 7^ 0. Plasmons in extrinsic 
semimetallic acGNR with the chemical potential in the lowest (highest) conduction (valence) band 
are found to have dispersion characteristics nearly identical to their intrinsic counterparts, with 
negligible differencs in dispersion arising from the slight differences in overlap integrals for the 
interband and intraband transitions. 



I. INTRODUCTION 

Plasmon propagationi"— in two-dimensional half-metallic 
graphene^^^ has been widely studied in both the Dirac 
(continuum) approximation^ and the orthogonal Pz- 
orbital tight binding (pzTB) approximatio n^ . In the 
latter, numerical integration of the Lindhard expression 
yields the dielectric function. 
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FIG. 1. (Color online) Ball and stick model of the hydrogen 
passivated acGNRS. The unit cell for this structure is high- 
lighted in gray, which consists of 16 carbon and 4 hydrogen 
atoms. X- and y- are the longitudinal and transverse direc- 
tions respectively. Atomic visualization is done by Hiickel- 

Armchair graphene nanoribbons (acGNRs)^"— of 
atomic width N of mod{N,3) = —1, where = 
2, 5, 8, 11, 14, 17, ... are also semimetallic within the con- 
tinuum and the pzTB approximation— i^i^^. One such 
acGNR of A^ = 8 with hydrogen passivation is shown in 
Fig. [TJ This semi-metallic behavior can be interpreted 
as a projection on the two-dimensional band structure of 
graphene, where the boundary conditions due to finite 
width of the ribbon crosses one of the two Dirac points, 
which are related by the time-reversal symmetry. 




FIG. 2. Band structure for acGNRs calculated using the 
pzTB model, (a) AT = 2, (b) AT = 5, (c) AT = 8, (d) AT = 11. 



Plasmons in the semi-metallic acGNRs with even 
atomic widths N = 2,8,14,... have been studied us- 
ing a continuum modeU^. Such acGNRs were shown to 
support plasmons that exhibited unique characteristics 
arising from the quasi one-dimensional character of the 
acGNRs. Further, we note that due to the nature of the 
boundary conditions required, acGNRs with odd atomic 
widths N = 5, 11, 17, ... are not analyzed using the con- 
tinuum model. Hence, a more detailed model is required 
to understand the plasmon dispersion characteristics of 
odd atomic width acGNRs. 

Plasmons in graphene microribbons have also been 
studied recently^^^^^. In one of these papers^^, the au- 
thors note that consideration of quantum confinement 
and microscopic structure of the nanoribbon becomes 
important for small ribbon widths and energies near the 
Dirac point. We address these issues in the present work. 

In this paper, we examine the plasmon dispersion in 
acGNRs^— by using a nearest-neighbor pzTB model 
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with a hopping parameter t = 2.beV ?^ These semimetal- 
hc acGNRs display hnear and symmetric dispersions for 
the low-lying conduction and valence bands within the 
pzTB theory. Although, when using a theory that in- 
cludes an extended basis set, the semimetallic acGNRs 
also develop a small bandgap on the order of a few meV 
and some nonlinearity around the band edge^^^^^, we ne- 
glect such effects in this work. For the cases where 
mod (A^, 3) = 0, 1, intrinsic acGNRs are semiconducting 
with appreciable bandgap and do not support plasmon 
propagation. 

New results presented here fall into four major ar- 
eas: 1) we present the first description of plasmons in 
odd-atomic-width semimetallic acGNR; 2) we show that 
the plasmon dispersion of all plasmons in semimetal- 
lic acGNR depends strongly on the localization of the 
Pz electronic wavefunctions; 3) we show that the con- 
duction/valence band overlap integral is zero only for 
k — k' =0, where k is the momentum of the initial state 
and k' is the momentum of the final state - implying that 
there is a decay path of non-zero probability for all plas- 
mons with g 7^ 0; and 4) where the plasmons from our 
pzTB model can be compared to the continuum results, 
we show quantitatively different dispersion relations. 



II. THEORETICAL MODEL 

The unit cell of the acGNR of atomic width N 
(acGNRA^) contains 2N atoms arranged in a honeycomb 
structure, as is shown in Fig. [1] for acGNRS. The 
unit vector is given as, a = dx^ where d = 3a cc and 
ttcc = 1.42 A is the carbon bond length. The pzTB 
Hamiltonian of this structure is a 2A' x 2A^ matrix con- 
taining only nearest-neighbor couplings. We transform 
the real-space Hamiltonain to the reciprocal space H{k) 
to calculate the eigen values Ei{k) and eigen functions 

c[^\k) for the eigen state i = 1, 2, 2A/', where i is 
also the band index and a represents the atomic loca- 
tion. For conduction bands, the band index ranges from 
i = 1 corresponding to the lowest lying conduction band 
and i = N corresponding to the highest lying conduction 
band. For valence bands, the index ranges from i = 1 for 
the highest lying valence band to i = N for the lowest 
lying valence band. 

The band structure for acGNRs with A/" = 2, 5, 8, 11 is 
shown in Fig. [21 the four narrowest acGNRs that exhibit 
linear conduction and valence bands with electron-hole 
symmetry. We note here that acGNRs with odd atomic 
widths A" = 5, 11, 17, ... are not considered in the contin- 
uum theory although they do exhibit the semi- metallic 
characteristics necessary to allow plasmon propagation. 

Next, we calculate the longitudinal dielectric function 
for acGNRs in the RPA, where setting the dielectric func- 
tion equal to zero gives the plasmon dispersion relation. 
In the RPA, the dielectric matrix for acGNRs can be 
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FIG. 3. Overlap integral \ {k'; l\k; 1)| for the lowest conduction 
band (1) and the highest valence band (1) as a function of k for 
acGNR2 (overlaps for the linear bands in other semimetallic 
acGNR exhibit similar behavior). In the top panel, we present 
a qualitative diagram of the nature of the overlap integral. In 
this panel, lighter shades indicate larger overlap and darker 
shades indicate smaller overlap. The shading of the upper 
left and lower right parts of the diagram (regions where the 
overlap integral ^ 1) has been adjusted to clearly show the 
small variation of the overlap from 0.97-1. The shading of 
the upper right and lower left parts of the diagram (regions 
where the overlap integral <^ 1) has been adjusted to clearly 
show the variation of the overlap from 0-0.15. In the bottom 
panel, a quantitative comparison of the overlap integral | {k' = 
37r/4:d\k)\ for the pzTB (solid curve) and continuum (dotted 
curve) models is shown. The pzTB curve depicts a section 
through the top panel at /c^d/vr = 0.75. It should be noted 
that the pzTB overlap is zero only for A; = A;^, in contrast to 
the result for the continuum model where the overlap is zero 
for all k, k' of the same sign. Thus, in contrast to results 
predicted from the continuum model, for nonzero values of 
q backscattering will be allowed. As a result, in the pzTB 
approximation there is a decay path for plasmons into free 
carriers for all ^ 7^ 0. 

written as^-, 

^ijmn{q,^) = SimSjn " Vijmn{q) Ilmn{q,Oj) (1) 

where Vijmniq) is the Coulomb matrix element in one 
dimension, UmniQ^ ^) is the polarizability of the acGNRs, 
and z, j, m, and n are the band indices. Non- trivial 
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FIG. 4. vj i j i(q) for acGNRA^ calculated using Eq. [12] with w = 1.08984 A and 0.90 A. Ribbon atomic widths for each panel 
are: (a) N — 2^ (b) N — '^^ (c) N — 8, and (d) N — 11. For comparison purposes, vj -^j ^{q) from Eq. 9 of Ref.— is shown as 
well for the two nanoribbons of even atomic width. 



solutions to the field equations requires, 

det [eijmn{q,^)] = (2) 

In the RPA for small the polarizability of an acGNR 
can be written as: 

Umn{q,^) = lim l^x 

(3) 

where m and n are band indices, ^5 = 2 is the spin de- 
generacy, Lx is the sample length, k is the momentum 
of the initial state, k' = k ^ q \s the momentum of the 
final state, and f{E) = !/[! + e^^"^)/^^^] is the Fermi- 
Dirac distribution function with chemical potential /i and 
Boltzmann constant /c^, where T is the temperature in 
K. h is the reduced Planck's constant and 7^ is a small 
number. We consider intrinsic acGNRs with the chemical 
potential /i = 0. The denominator of the polarizability 
diverges for the energy difference between initial and final 
states equal to the plasmon energy hu. 



The unit cell for the nanoribbon is two-dimensional, 
however because d <^ q~^ for small q^ we write the ma- 
trix element of the Fourier transform of the Coulomb 
potential as: 

Vijmn {q) = ^ JJJ j (l)i{k;y,x)(l)j{k^q;y, x) 
xKo{q\y-y'\) 

X ^^(/c + q; y\ x) (pnik] y\ x) dy dy' dx dx 

(4) 

where the integral is taken over the entire two- 
dimensional unit cell, and where KQ{q\y — y'\) is the 
modified Bessel function of the second kind of order 0. 
For acGNRs, the above integral is independent of k^. 
(j)i{k]y^x) are eigenfunctions of the i-th band. In the 
pzTB, we expand the eigenfunctions by writing the basis 
as a set of ^-functions: 

0* (/c; y, x) (t)j {k + x) = ^ 4"^* (^) 4"^ + ^) 

X ^{y - yoc)^{x - Xa) (5) 

where the sum over a runs over all atoms in the unit 
cell, (ya^Xa) is the location of the a'th atomic site, and 
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FIG. 5. Various off-diagonal v{q) matrix elements for acGNR5 
calculated using Eq. [12] with w — 1.08984 A. The magnitude 
of the matrix elements is shown in the top panel, and the 
phase of the matrix elements is shown in the bottom panel. 



the c-"^(/c)'s are the complex eigenvector coefficients at 
each atomic site. Substituting this expansion into Eqn. 
m it can easily be shown that VijmniQ) separates into 
two factors Y{q) X{q) where Y{q) is the integral over the 
transverse {y) coordinate and X{q) is the integral over the 
longitudinal {x) coordinate. In the following, we solve for 
each of these two factors separately, and then combine 
them to obtain VijmniQ)- 

Substituting Eqn. [5] into Eqn. |4l the transverse factor 
becomes: 



xKo{q\ya-y0\) c'i>{k + q)ci^\k) (6) 



The ya = y(3 terms in Eqn. [6] diverge due to the di- 
vergence of the Bessel function for small arguments. We 
generalize the basis set for these terms, to take advantage 
of the convergence of the screening integral for basis sets 
of non-zero width, and obtain for the y^ = y^ terms of 
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f^°'\k;q) j j veci{y)Ko{q\y-y'\)veci{y')dydy' (7) 



with: 



=cf)*(fc)c(% + g)c(;^)*(fc + g)c(r)(fc) (8) 

and: 

«(y + f)-«(y-f) 



rect(y) 



(9) 



where u{y) is the Heaviside theta function and w repre- 
sents the transverse {y) spatial extent of the Pz orbital. 
The longitudinal factor X{q) is written: 



d d 



X{q) = J j 5{x-Xa)5{x' -x^)dxdx' (10) 







We set X = x' io compute the one-dimensional screening 
result for VijmniQ) in the d <C q~^ limit. In this limit, 
X{q) diverges; however by generalizing the (5- functions as 
rect(x) functions defined by Eqn. [9l we obtain: 



X{q) 



Xa Xp 



(11) 



Finally, combining Eqns. [6land[TT| we find: 
2e^ 3a^, 



eo w /-^ 

xK,{q\y^-yp\) c'i^ik + q) ci^\k) (12) 



We note that the magnitude of Vijmn{Q) (and therefore, 
the plasmon dispersion) is quite sensitive to the value 
of w for two reasons: 1) the factor 3acc/w in Vijmn{Q) 
scales the magnitude directly; and 2) the dependence of 
the ya = yf3 terms of VijmniQ) on w through Eqn. [71 
scales the shape of the Bessel integral for small q. As w 
decreases, VijmniQ) increases and as a result, the phase 
velocity of the pzTB plasmon also increases. 

Two-Band Model: For transitions other than those 
that involve the highest lying (linear) valence band and 
lowest lying (linear) conduction band, the Coulomb ma- 
trix elements for non-symmetric bands are small and the 
energy differences are large relative to those for the lin- 
ear bands as discussed in the next section. As a re- 
sult, contributions to the dielectric function from bands 
other than the linear bands may be neglected, and a 2- 
band dielectric function that includes contributions from 
only the hnear bands represents a good approximation 
to the complete 2A/'-band dielectric function. Therefore, 
we are interested in plasmons for q ^ 0, where contri- 
butions to the polarizability are dominant for transitions 
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FIG. 6. Plasmon dispersion curves for acGNRA/" calculated using Eq. [13] with w — 1.08984 A and 0.90 A. Ribbon atomic 
widths for each panel are: (a) iV = 2, (b) N — (<z) N — 8, and (d) N — W. For comparison purposes, the 2-band plasmon 
dispersion calculated at T = using Eqns. 14 and 15 of Ref.^^ is shown as well for the two nanoribbons of even atomic width. 
The free carrier dispersion is shown as a dotted line. 
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FIG. 7. Plot of the plasmon wavenumber q normalized by 
the free space wavenumber as a function of free space 
wavelength Ao for wavefunction locahzation parameter 
w 1.08984 A. 




FIG. 8. The group velocity Vg of acGNRAT (N = 2, 5, 8, 11) 
with ^ ^ is plotted for values of the wavefunction local- 
ization parameter near w — \K. 



between the highest lying (linear) valence band (m = 1) 
and the lowest lying (linear) conduction band (n = 1). 
The polarizabilities arising from transitions from lower 
valence bands to higher conduction bands are negligible 
and hence are neglected in our calculation of the plasmon 
dispersion relations. 



In intrinsic acGNRs at T = 0, the self-polarizabilities 
of the linear bands are given as, ^^^{(i^uj) = nii(g,cc;) = 
0. Further, symmetries in the acGNRs require^^^^^ that 
the Coulomb matrix elements '^i i y i (^) = '^i i i t(^) = 
'^1 T 1 t(^) ~ '^iTTi(^)- This result gives the dispersion 
relation of the collective (plasmon) state in the 2-band 
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FIG. 9. Schematic illustration of the difference between the 
intrinsic and extrinsic plasmon excitations. The linear con- 
duction and valence bands are shown as the diagonal bold 
lines. Intrinsic transitions are interband transitions that oc- 
cur across the intrinsic chemical potential /x = at the Dirac 
point from a filled (empty) state in the valence (conduction) 
band to an empty (filled) state in the conduction (valence) 
band. Intrinsic transitions to a higher (lower) energy state 
contribute to the IIj^ (^ii) polarizability term. Extrinsic 
transitions are intraband transitions that occur across the ex- 
trinsic chemical potential /i = /Hext from a filled (empty) state 
in the conduction band to an empty (filled) state in the con- 
duction band. All extrinsic transitions contribute to the IIn 
polarizability term. Negligible differences in the magnitude of 
the wavefunction overlaps (the overlaps are nearly 1 in both 
cases) between the interband and intraband transitions ac- 
count for the small group velocity differences between plas- 
mons. Otherwise the plasmon dispersion is identical between 
the two cases. 

approximation by simplifying Eq. 2 as follows, 

l-^T,i,T,i(^)[nTifec^) + n,y(g,c^)] =0 (13) 

III. DISCUSSION OF RESULTS 

When calculating the dielectric function of acGNRs for 
g 0, the overlap integral {k']j\k]i) plays an important 
role, where k' = k -\- k and k' are the initial and final 
carrier wave numbers, q is the plasmon wave number. We 
note that the overlap integral is identically zero for all 
but the bands symmetric about /i = 0. For example, the 
wavefunction overlap between the highest lying (linear) 
valence and lowest lying (linear) conduction bands for 
acGNR2 is shown in Fig. [3l The overlap integral is about 
unity for states k ^ —d/ir and k^ ^ d/n on opposite sides 
of the Brillioun zone (BZ). The value decreases linearly 
from unity for states immediately adjacent to the zone 
center, to about 0.97 for states k ^ ±7r/d and k^ ^ ^ir/d. 

For overlaps between states on the same side of the 
BZ the overlap is identically zero for k = k' and varies 



linearly up to about 0.15 for states near opposite edges of 
the half-BZ. This latter result indicates that optical (free- 
carrier) transitions with small, non-zero q are allowed in 
the pzTB formalism, and can be expected to have an 
impact on the Landau damping of plasmons with q 0. 

For higher lying bands, transitions between states on 
the same side of the BZ are qualitatively similar, with 
k = k' transitions forbidden, and k ^ k' transitions al- 
lowed. For plasmon transitions through the zone cen- 
ter, the overlap is more detailed. However, because the 
energy difference between states in the higher bands is 
much larger than that of the highest valence band to 
lowest conduction band transition, such transitions do 
not contribute significantly to the dielectric function and 
are neglected in this paper. 

As noted above, VijmniQ) is quite sensitive to the lo- 
calization of the Pz wavefunctions at each atomic site. 
Calculations using extended Hiickel theory^^ suggest the 
localization is approximately w = 1 A. One interesting 
choice for w is to require that the area of the generalized 
atomic wavefunction S{y) S{x) Tect{y) rect(x) be equal 
to the area of the circular disk surrounding each atomic 
site with radius r = >/3acc/2 (disks are tangent to the 
bisector line between longitudinal rows of atoms in the 
acGNR). This results in a localization parameter w = 
1.08984 A and gives remarkable agreement between the 
plasmon dispersion relation calculated using the pzTB 
and continuum models as discussed below. In Fig. |4]we 
plot the pzTB Coulomb matrix elements vj ^j-^ (q) of Eq. 
[12] for AT = 2,5,8,11 using w = 1.08984 A' and 0.90 A. 
For reference, we also plot '^i i y i(^) for the continuum 
model (Eq. 9 of Ref.^^) for '7v '= 2 and 8. This ma- 
trix element describes the coupling between the highest 
lying (linear) valence band and the lowest lying (linear) 
conduction band. In the small q limit, the matrix ele- 
ments calculated using the continuum model agree with 
the pzTB results for w = 1.08984 A. Vijmn{Q) is seen 
to increase as w decreases. Dispersion of the pzTB ma- 
trix element is larger for small q, and as a result the ra- 
tio of the pzTB and continuum matrix elements depends 
strongly on the acGNR width as q increases. 

For 2M-band dielectric functions with M > 1, 
Coulomb matrix elements that describe coupling between 
higher lying and lower lying bands also appear in the 
expression for the dielectric function. We plot the mag- 
nitude and phase of several of these additional matrix 
elements in Fig. [5l In the small q limit, the magni- 
tude of each of these matrix elements goes to zero , and 
their magnitude remains small relative to the magnitude 
of vj ^ even as q increases. Such a result confirms the 
validity of using a 2-band approximation to the dielectric 
function when calculating the plasmon dispersion in the 
small q limit. 

Next, we highlight several differences between the con- 
tinuum model and the pzTB model. Perhaps the most 
significant of those differences is the existence of plas- 
mons for acGNRs with odd atomic width A" = 5, 11, 17, ... 
Results for the 2-band pzTB plasmon dispersion of 
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acGNRs with = 2,5,8 and 11 calculated using Eq. [T]in 
the 2-band approximation are shown in Fig. [6l Because 
the plasmon dispersion is quite sensitive to the particular 
localization parameter chosen, we show the dispersion for 
two values of w = 1.08984 A and 0.90 A. For comparison, 
we also plot the plasmon dispersion calculated using the 
2-band continuum model. Here we emphasize that the 
continuum 2-band dispersion relation does not apply for 
A/' = 5, 11. This results from the fact that the continuum 
analysis applies only to acGNRs with even atomic width, 
due to the boundary conditions on the electronic wave- 
functions applied in the continuum model. Although the 
agreement between the phase velocities of the plasmons 
from the two different models for w = 1.08984 A is re- 
markable, it is difficult to attach physical signficance to 
this result, since value of w is not known with accuracy 
better than w ^ 1 A. 

In Fig. [7] we plot the free-space-normalized plas- 
mon wavenumber q vs. free-space wavelength Aq for 
a value of the pz wavefunction localization parameter 
w = 1.08984 A. This fi gure illustrates the wavelength 
dependence of the plasmon phase velocity. 

Several differences are noted when comparing the 
pzTB and continuum results. First, the group velocity of 
the plasmons for ^ is smaller for the continuum model 
than for the pzTB result. Second, the magnitude of the 
group velocity dispersion of the pzTB plasmons is greater 
than that of the continuum model plasmons for small q. 
These differences arise primarily from the greater disper- 
sion of the Coulomb matrix element vj -^j {q) at small 
q for the pzTB model than for its continuum counter- 
part. We note that the pzTB plasmon dispersion de- 
pends strongly on the wavefunction localization through 
the parameter and as a result for wavefunctions that 
are more strongly localized, the pzTB plasmons have a 
larger phase and group velocity than their continuum 
counterparts as well. This dependence of group velocity 
on for g ^ is shown in in Fig. [51 

Extrinsic acGNR: Here we consider the plasmon dis- 
persion of electron-doped extrinsic semimetallic acGNR 
where the chemical potential fi = fiext is in the i = 1 
conduction band. At T = 0, this corresponds to filled 
electronic states in the conduction band between —kF< 
k < kp where kp = /^ext/^VF- For this case, the plasmon 
dispersion arises from intraband transitions in the i = 1 
conduction band. Due to symmetries inherent in the 
acGNR, = vj^ij^iiq)' Further, Uii{q,uj) ^ 0, 

and for plasmons with q < kp^ Hy-^ and H-^y ^ 0. With 
these simplifications, the plasmon dispersion of Eqn. [2] 
reduces to: 



l-^T,i,T,i(^)niifec^) = (14) 



with: 

Hii((7,a;) = 

(15) 

and: 

A(fc,^) = |^fc+,,i|-|^fci| (16) 

For the intraband transition, the overlap integral {k + 
q]l\k]l) ~ 1 for k -\- q and k of the same sign, and is 
nearly otherwise. As a result, the extrinsic plasmon 
dispersion of Eqn. [ll]is identical to the intrinsic plasmon 
dispersion of Eqn. [13] except for small differences in the 
value of the overlap integral used in the calculation of the 
relevant polarizability. For the extrinsic case, the tran- 
sitions contributing to the plasmon dispersion are across 
the chemical potential /i (on the same side of the BZ), 
whereas for the intrinsic case, the transitions contribut- 
ing are across the Dirac point. A schematic illustration 
of these two types of processes is shown in Fig. [9l 

Extrinsic plasmons with q > kp include both intraband 
and interband contributions to the dielectric function. 
In this case the dispersion relations are identical as well 
due to the similarity of the polarizabilities, except for 
the small differences in overlap between the two types of 
transitions. 

Arguments similar to the above also show that ex- 
trinsic plasmons in semimetallic acGNR with a negative 
chemical potential will have nearly identical dispersion 
to the intrinsic plasmons. 

IV. EXPERIMENTAL CONSIDERATIONS 

In this section, we discuss briefly experimental con- 
siderations for demonstrating the propagation of acGNR 
plasmons in the THz regime. First, we address possi- 
ble coupling mechanisms for exciting the plasmon modes. 
Efficient coupling into acGNR plasmon modes is difficult 
to achieve due to momentum-conservation issues - the 
nanoribbon plasmon momentum is much larger than the 
photon momentum of the excitation source (see Fig. [7]). 
There are two primary types of excitation mechanisms 
that one may consider: 1) point-source excitation and 2) 
distributed excitation. Much work has gone into the de- 
velopment of point-source THz sources, motivated by the 
initial work of Katzenellenbogen and Grischkowsky^^. In 
principle, a point-source of this type, coupled directly 
to an acGNR should be able to efficiently excite a THz 
plasmon. More recently, experimental work on plasmon 
excitation using distributed light-plasmon coupling^^ has 
been successful at exciting plasmons in microribbon ar- 
rays. This experimental approach was used to study plas- 
mon resonances in the microribbons. 
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A second consideration for a successful experiment is 
the substrate on which the acGNR wih be mounted. If 
the substrate has a different dielectric constant than the 
overlayer in which the acGNR is embedded, the presence 
of the substrate dielectric material must be taken into 
account. An effective dielectric constant may be defined 
that is a weighted spatial average over the dielectric re- 
gions surrounding the acGNR. It is this effective dielec- 
tric constant that will appear in Eqn. |4] for %mn(^)- 
Also, the substrate dielectric material may have some 
losses associated with it due to non-zero conductivity a 
that must be taken into account. 



The issue of plasmon losses is also important in a 
carrier-scattering context. In a recent experimental 
paper^^, the carrier dephasing time in graphene has 
been measured for excitation energies below the optical 
phonon energy hwopt = 200 meV. Measured values were 
nominally T2 = 25 ps. Using this value for r2, assuming 
a plasmon velocity of Vg = 4.5 x 10^ m/s (see Fig. [8]) and 
a plasmon wavevector q = 0.005 A~^, it can be shown 
that the propagation distance for the plasmon may be as 
large as ~ 1000 wavelengths. 



V. SUMMARY 

In this paper, we have calculated the plasmon disper- 
sion for semimetallic acGNRs using a pzTB model. Con- 
tributions to the dielectric function from the linear con- 
duction and valence bands are included. Coulomb matrix 
elements for other transitions are small and go to zero as 
g ^ 0, and the energy differences due to higher lying 
bands are significantly larger than from the linear bands. 
As a result, contributions to the plasmon dispersion from 
these higher transitions are neglected. The results show 
several differences in the plasmon dispersion from that 
obtained using a continuum model, and include a group 
velocity for small q that depends strongly on the Pz wave- 
function localization, and a larger group velocity disper- 
sion. Further, and perhaps most significantly, the pzTB 
model predicts the existence of plasmons in acGNRs with 
odd atomic width N that were not identified from the 
continuum model. Finally, the overlap integrals for tran- 
sitions with small, non-zero q in acGNRs are shown to 
be small but non-zero. This result suggests that such 
plasmons should exhibit Landau damping, although the 
damping rate is expected to be small. 
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